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Abstract 

We apply the /ip-version of the boundary element method (BEM) for the numerical so- 
lution of the electric field integral equation (EFIE) on a Lipschitz polyhedral surface T. The 
underlying meshes are supposed to be quasi- uniform triangulations of T, and the approx- 
imations are based on either Raviart-Thomas or Brezzi-Douglas-Marini families of surface 
elements. Non-smoothness of T leads to singularities in the solution of the EFIE, severely 
affecting convergence rates of the BEM. However, the singular behaviour of the solution can 
be explicitly specified using a finite set of power functions (vertex-, edge-, and vertex-edge 
singularities). In this paper we use this fact to perform an a priori error analysis of the hp- 
BEM on quasi-uniform meshes. We prove precise error estimates in terms of the polynomial 
degree p, the mesh size h, and the singularity exponents. 
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1 Introduction 

In this paper we study numerical approximations of the electric field integral equation (EFIE) 
on a surface T = 90, where f2 C M 3 is a Lipschitz polyhedron. The EFIE is a boundary integral 
formulation of a boundary value problem for the time-harmonic Maxwell equations in the domain 
exterior to O. It models the scattering of electromagnetic waves at a perfectly conducting body. 

For the numerical solution of the EFIE we use the Galerkin boundary element method 
(BEM). It employs H(divr)-conforming families of surface elements, namely Raviart-Thomas 
(RT) and Brezzi-Douglas-Marini (BDM) elements, to discretise the variational formulation of 
the EFIE (called Rumsey's principle). This approach is referred to as the natural BEM for 
the EFIE. As in finite element methods, the convergence may be achieved either by keeping 
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polynomial degrees fixed and refining the mesh (h- version), or by fixing the mesh and increasing 
polynomial degrees (p- version), or by simultaneous /i-refinement and p-enrichment (hp- version). 

The Galerkin BEM is widely used in the engineering practice for the simulation of electro- 
magnetic scattering. Moreover, it had been used long before a rigorous theoretical analysis of 
the method became available. Error analyses of different boundary element schemes for the eddy 
current problem on smooth obstacles are given by MacCamy and Stephan in [2i\ 125] . Despite 
these early results the BEM-analysis for non-smooth obstacles started much later. In fact, the 
convergence and a priori error analysis of the /i-BEM for the EFIE on piecewise smooth (open 
or closed) surfaces has been developed within the last decade (see [131 [22l [HJ [101 EE] ) , and the 
corresponding results for high-order methods (p- and /ip-BEM) are very recent (see [6j SI E] ) - In 
particular, [6] and [7] establish quasi-optimal convergence of the natural /ip-BEM on meshes of 
shape-regular elements. An essential ingredient for the proofs in these papers are the projection- 
based interpolation operators developed by Demkowicz and co-authors [El [20J, [18] . These oper- 
ators also facilitate a priori error analysis, where one needs hp- approximation theory in specific 
trace spaces. In [4j we prove an a priori error estimate of the /ip-BEM on quasi- uniform meshes 
of affine elements under the assumption that the regularity of the solution to the EFIE is given 
in Sobolev spaces on V. The latter result states that the method converges with the same rate 
r + i in both h and p^ 1 (here, r denotes the Sobolev regularity order, and p is assumed to 
be large enough). However, it is evident from the numerical results reported in [23J that the 
p-BEM converges faster than the /i-BEM for the EFIE. This is similar to what was observed 
and proved for the p-BEM and the /i-BEM applied to elliptic problems in three-dimensions (cf. 

[221 EH m [5]). 

With this paper we fill a gap in the theory of the BEM for the EFIE on polyhedral surfaces 
by proving a precise error estimate for the /ip-BEM on quasi- uniform meshes. Similarly to 
the elliptic case we make use of explicit expressions for singularities in electromagnetics fields 
(these expressions are available from [171 [6]). The established error estimate shows that the 
convergence rates in h and p depend on the strongest singularity exponent, and that the p-BEM 
converges twice as fast as the /i-BEM on quasi-uniform meshes. This extends the results of [6], 
where the analysis was restricted to the p-BEM on a plane open surface. 

It is now well known that appropriate decompositions of vector fields (on both the con- 
tinuous and discrete level) are critical for the analysis of electromagnetic problems and their 
approximations. In the case of the EFIE the main idea is to isolate the kernel of the divr- 
operator such that the complementary field possesses an enhanced smoothness (cf . [101 El Ej ) . 
The corresponding decompositions of singular vector fields greatly facilitate our analysis in this 
paper as well. In particular, the complementary vector fields are singular vector functions be- 
longing to H 1 / 2 (r). Then, in the case of the BDM-based BEM, these vector functions can be 
approximated component-wise by continuous piecewise polynomials, and the desired hp-error 
estimates are derived by using the corresponding results in [3] for scalar singularities (belonging 
to H l l 2 (Y)). However, this simple approach does not work for the RT-based BEM because the 
dimension of the underlying RT-space on the reference element is smaller than the dimension of 
the BDM-counterpart. That is why, the results of [3j cannot be applied directly, and additional 
technical arguments are needed (see the proof of Lemma 14.21 in Section f5.2j) . 



2 



The rest of the paper is organised as follows. In the next section we formulate the EFIE 
in its variational form and recall the typical structure of the solution to this model problem. 
In Section [3] we introduce the hp-version of the BEM for the EFIE and formulate the main 
result of the paper (Theorem 13. ip stating convergence rates of the method. This result follows 
from the general approximation theorem (Theorem 14. 1 j) established in Section 2J The proof of 
Theorem 14 . 1 1 relies . in particular, on two technical lemmas which are proved in Section [5j 

Throughout the paper, C denotes a generic positive constant independent of h and p. 

2 Formulation of the problem 

Let r be a Lipschitz polyhedral surface in M 3 . Throughout the paper we will use exactly the 
same notation as in [3] for all involved differential and boundary integral operators as well as 
for Sobolev spaces of scalar functions and tangential vector fields on T (all essential definitions 
are given in [U Section 3.1]). In particular, we use boldface symbols for vector fields, and the 
spaces (or sets) of vector fields are denoted in boldface as well. 

For a given wave number n > 0, we denote by \& K (resp., \I/ K ) the scalar (resp., the vectorial) 
single layer boundary integral operator on T for the Helmholtz operator — A — k 2 (see |15l 
Section 5]). The variational formulation for the EFIE will be posed in the following Hilbert 
space of tangential vector fields on T: 

X = H-V^divr,! 1 ) : = { u G *V /2 (r); div r u G H-V 2 (T)}, 

which is the trace space of H(curl, Q), where fl C M 3 is a Lipschitz polyhedron such that F = dQ 
(we refer to [THESE] for the definition and properties of this and other trace spaces on T). 

Let X' be the dual space of X (with duality pairing extending the L 2 (r)-inner product for 
tangential vector fields). Then, for a given source functional f £ X' the variational formulation 
of the EFIE reads as: find a complex tangential field u € X such that 

o(u,v) := (^ K div r u,div r v) -k 2 (* k u,v) = (f,v) Vv € X. (2.1) 

Here, the brackets (■,•) denote dualities associated with H l l 2 (T) and Hy^ 2 (r). To ensure the 

uniqueness of the solution to (12. ip we always assume that k 2 is not an electrical eigenvalue of 
the interior problem. 

Let us recall the typical structure of the solution u to problem (|2.ip , provided that the source 
functional f is sufficiently smooth (we note that this regularity assumption is satisfied for the 
electromagnetic scattering with plane incident wave). We use the results of [6, Appendix A]. 
These results were derived from the regularity theory in [17J for the boundary value problems 
for Maxwell's equations in 3D by making use of trace arguments and some technical calculations 
in the spirit of [29] and |28| . 

Let V = {v} and E = {e} denote the sets of vertices and edges of T, respectively. For v € V, 
let E(v) denote the set of edges with v as an end point. Then the solution u of (12. ip can be 
written as 

U = U rcg + U sing , (2.2) 
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where 

u reg G X fe := {u G H fc _(r); div r u G H k _(T)} with k > (2.3) 

(the spaces H*L(r) and H k _(T) are defined in a piecewise fashion by localisation to each face of 
T, and the space X fc is equipped with its graph norm || • || X fc, see [1]), 

usin g = E ue + E uU + E E ^ ( 2 - 4 ) 

and u e , u v , and u ei; are the edge, vertex, and edge- vertex singularities, respectively. 

In order to write explicit expressions for the above singularities, let us fix a vertex v G V 
and an edge e G E(v). Then, on each face T ev C V such that e C dT ev we will use local polar 
and Cartesian coordinate systems (r v ,9 v ) and (x e i,x e 2), both with the origin at v, such that 
e = {(x e i, x e 2)\ x e 2 = 0, x e i > 0} and for a sufficiently small neighbourhood B T of v there holds 
T ei; n B T C {(r\,,6^); < 9 V < lo v }. Here, w„ denotes the interior angle (on T ev ) between the 
edges meeting at v. For simplicity of notation we write out here only the leading singularities 
in u e , u v , and u ev on the face T ev , thus omitting the corresponding terms of higher regularity 
(see [BJ Appendix A] for complete expansions). 

For the edge singularities u e one has 

u e = curlr- (xll |logx e2 | s i b\ (x el ) Xi(^el) xK 3 ^)) + x£ |logx e2 | s 2 bf(x e i) Xii^d) xK 3 ^), 

(2-5) 

where curlr-™ = (d/dx e 2, —d/dx e i) is the tangential vector curl operator curlr restricted to 
the face T ev (cf. [HIE]), jf, 7| > ^, and sf, s| > are integers. Here, x\-, X2 are cut-off 
functions with x\ = 1 m a certain distance to the end points of e and xf = in a neighbourhood 
of these vertices. Moreover, x% = 1 f° r < £ e 2 < <5 e and X2 = for x e2 > 2<5 e with some 
5 e G (0, |). The functions 6fxf G H mi (e) and b|xf G H m2 (e) for mi and 777-2 as large as 
required. 

The vertex singularities u v have the form 

= curlr«(r^ |logr„|tf x v (r v )xWv)) + r$ \]ogr v f> X v (rv)xl(0v), (2.6) 

where A^, A2 > — \ are real numbers, qf, q^ > are integers, x^ is a C°° cut-off function with 
X v = 1 for < r v < t v and x" = for r„ > 2r^ with some t„ G (0, ^). The functions Xi\ X2 are 
such that xf G .fZ"* 1 (0, w„), X2 e H* 2 (0,Wt,) for ij, £2 as large as required. 
For the combined edge-vertex singularity u 6 " one has 

where 

uf = curlre.(x e A r 7f 4|log a ; el |^|log a ; e 2| /32 x t, (^)x ei '(^)) 
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+ x2 ^xll\\ogx el \^\\ogx e ^ X v {r v )x ev {Ov) ( " J (2.7) 

and 

uf = curlpe^x^ | \ogx e2 \ s ° b%{x e i,x e2 ) xl^ei)) + %e2 I logx e2 | s 2 b|(x e i,x e2 ) xK^)- (2.8) 

Here, A", 7?, sf (i = 1, 2), x"i an d X% are as above, ^ > (k = 1 . . . , 4) are integers, f3± + /3 2 = 
sf + ?i , Pz + & = s| + q v 2 with gf being as in f^H]) . x CT is a C°° cut-off function with x e " = 1 
for < 9 V < (3 V and \ ev = f° r §Aj < < Wd for some /3„ € (0, min{^/2, 7r/8}]. The functions 
6| and b|, when extended by zero onto IR 2+ := {(x e i,x e 2); x e2 > 0}, lie in H mi (SR 2+ ) and 
H m2 (IR 2+ ), respectively, with mi, m 2 as large as required. Finally, the supports of uf u and 
are subsets of the sector S ev = {(r v ,0 v ); 0<r v < 2t v , 0<0 v < §/3„}. 

Remark 2.1 (i) The exponents 7 4 e (i = l,2) of the edge and vertex-edge singularities in (|2.5|) , 
(12771) . (USD sate/?/ 7 4 e > \. However, for our approximation analysis below it suffices to require 
that 7? > (i = 1,2). Note that 7 4 e > and A^ > — \ (i = 1,2) are the minimum requirements 
to guarantee u € X. 

(ii) As mentioned above, the terms of higher regularity (i.e., with greater singularity expo- 
nents) are omitted in (|2.5p - (|2.8p . These terms are necessary to obtain the regular part u reg £ X fc 
of decomposition (|2.2h as smooth as required. This can be done by considering sufficiently many 
(omitted) singularity terms of each type. 

Remark 2.2 (i) By (j2.5|) ()2.8j) we conclude that any singular vector field u s in (|2.4p (s = e, v, 

or ev ) can be written as 

u s = curl r w s + v s = curl r w s + (vf, uf ) (2.9) 

with corresponding (scalar) singular functions w s , v\, t>| being defined on the whole surface T. 
Note that these scalar functions are H 1 / 2 -regular on T, i.e., 

w'eH^iT), v s = (v{, V2) € H 1/2 (r), s = e, v, ev. 

(ii) It is important to observe that the functions w s , vf , v 2 (s = e, v , or ev) in (12.9P are typical 
scalar singularities inherent to solutions of the boundary integral equations with hypersingular 
operator for the Laplacian on V and with possibly singular right-hand side. Continuous piecewise 
polynomial approximations of these scalar singularities in fractional- order Sobolev spaces were 
analysed in [3j [l] and will be used to prove the main result of the present paper. 

3 The hp- version of the BEM and the main result 

For the approximate solution of f)2. 1 1) we apply the /ip-version of the BEM on quasi-uniform 
triangulations of T. Our BEM is based on Galerkin discretisations with an appropriate family of 
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H(divr, r)-conforming surface elements (here, H(divr,T) := {u £ L 2 (r); divr u € L 2 (T)}). In 
what follows, h > and p > 1 will always specify the mesh parameter and a polynomial degree, 
respectively. For any C M. n we will denote pn = sup{diam(i3); B is a ball in Vt\. 

Let T = {A^} be a family of meshes A/, = -{Tj; j = 1, . . . , J}. Each mesh is a partition of 
r into triangular elements Tj such that f = Uj =1 Tj, and the intersection of any two elements 
Tj> ffc (j ^ A;) is either a common vertex, an entire side, or empty. 

In the following we always identify a face of the polyhedron T with a subdomain of M. 2 . 
We denote hj = diam(Fj) for any Tj € A/j. Furthermore, any element Tj is the image of the 
reference triangle K = {(£i, £2); < £1 < 1, < £2 < 1 — £1} under an affine mapping Tj, more 
precisely 

tj = Tj(K), x 7' i£; £ • bj, 

where Bj € M 2x2 , bj € M 2 , x = (xi,x 2 ) € f j, and $ = (£1^2) € Then, the Jacobian matrix 
of Tj is Bj £ l 2x2 , and its determinant Jj := det (-Bj) satisfies the relation \Jj\ ~ h 2 . 

We consider a family T of quasi-uniform shape-regular meshes A^ on T in the sense that 

there exist positive constants a%, cr 2 independent of h = maxL such that for any Tj G A^ and 

i 

arbitrary A/j £ T there holds 

hj<aipr 3 , h<a 2 hj. (3.1) 

Whereas the mapping Tj introduced above is used to associate scalar functions defined on the 
real element Tj and on the reference triangle K, the Piola transformation is used to transform 
vector- valued functions between K and Tj\ 

v = Mj(v) = j-Bji? o Tf\ v = MjV) = J i^7 lv T i- ( 3 - 2 ) 

Let us introduce the needed polynomial sets. By V P (K) we denote the set of polynomials of 
total degree < p on the reference triangle K. We will use two families of H(divr, r)-conforming 
surface elements: the Raviart-Thomas (RT) and Brezzi-Douglas-Marini (BDM) elements. The 
corresponding spaces of degree p > 1 on the reference triangle K will be denoted as follows (see, 
e.g., UM)- 

Vf{K) = (V^iK)) 2 ZV p -i{K); 
V^ U {K) = (Vp(K)) 2 . 

We will use the unified notation "P V {K) which refers to either the RT- or BDM-space on K for 
p > 1. Accordingly, all results in this paper are formulated in a unified way and are valid for 
both the RT- and BDM-based boundary element spaces defined on the triangulation of T. Note, 
however, that in some cases we will need to provide arguments separately for each type of these 
boundary elements. 

Using transformations (|3.2h . we set 

X hp :={vGH(div r ,r); Mj\v\ T .) £V P {K), j = 1, . . . , J}. (3.3) 
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Note that only one type of surface elements (i.e., either the RT- or BDM-elements) is used in 
(|3.3p for all triangles Tj. We will denote by N = N(h,p) the dimension of the discrete space 
X/^. One has N ~ h~ 2 for fixed p and iV ~ p 2 for fixed h. 

The hp-version of the Galerkin BEM for the EFIE reads as: Find G X/j p such that 

a(u hp ,v) = (f,v) Vv G X. hp . (3.4) 



Due to the infinite-dimensional kernel of the divr-operator, the bilinear form a(-, •) in (|2.ip 
is not X-coercive, and, hence, the unique solvability of (13.41) cannot be proved by standard 
arguments. However, the refined analysis in [7| shows that the unique BEM-solution u^ p G 
X/jp does exist, and it converges quasi-optimally to the exact solution u € X of the EFIE as 
N(h,p) — > oo. This result is formulated in the following proposition. 

Proposition 3.1 Theorem 1.2] There exists No > 1 such that for any f 6 X' and for 

arbitrary mesh-degree combination satisfying N(h,p) > Nq the discrete problem (j3.4j) is uniquely 
solvable and the hp-version of the Galerkin BEM converges quasi-optimally, i.e., 

||u - u^pHx < Cinf{||u - v|| x ; v G X fep }. (3.5) 

Here, u G X is the solution of (|2.ip . G X/ ip is the solution of ()3.4|) . || • ||x denotes the norm 
in X, and C > is a constant independent of h and p. 

The next theorem is the main result of this paper. It states convergence rates of the np-BEM 
with quasi-uniform meshes for the EFIE. These convergence rates (in both the mesh parameter 
h and the polynomial degree p) are given explicitly in terms of the singularity exponents of the 
vector fields in {22HZ2]) • 

Theorem 3.1 Let u G X and Uh p G X/j P be the solutions of (12. ID and (13. 4p . respectively. 
We assume that the source functional f in f)2. 1 1) is sufficiently smooth such that representation 
(I2.2j) - (l2.8p holds for the solution u of (|2.ip . Let vq G V, eo G E(vo) be a vertex-edge pair such 
that 

min{\? + 1/2, A* + 1/2, 7l e ° , 72°} = min min{A? + 1/2, A?J + l/2, 7l e ,7 2 e } 

u?i£/i A^ and 7 4 e = 1,2,) 6eing as in (I2.5p - (l2.8p . T/ien /or any /i > and /or every p > 
min{Af ,A^, 7l e ° - ±, 7 |° - |} i/iere ZioWs 

/ , x min{ A? +1/2, A^ +l/2, 7l e ° , 7 |0 } 

||u-u,p||x<C^J + ' (3 ' 6) 



where 



{ max {of + sT + h 12 + *2° + 3> # A *° = 7f° - | /or i = 1, 2 
max {<° + S f + i, + s?} if Af = 7l eo - \, X v 2 ° + 72 e 



5 ^; 

eo 



max {gf + s * <g° + s *> + 1} if A?° + 7 f - I, A^ = 72 e ° - ± 

„ max {q^° + q^ + s 2°} otherwise 



(3.7) 



2' 



7 



with the numbers and qj° (i = 1,2) given in (|2.5[) and (|2.6p . respectively, and 

\ ifp = min{Af , A^, 7l eo - l/2,7 2 e ° - 1/2}, 
otherwise. 



(3.8) 



Ifl<p< min {Ai°, A2 , 7i° - |, 72° - \}, then for any h > t/iere /iotas 

||u-u^||x<C/i p+1 / 2 . (3.9) 



Proof. Considering enough singularity terms in representation (|2.4p the function u reg in (|2.3p 
is as regular as needed. Then, due to the quasi-optimal convergence (|3.5h of the /ip-BEM with 
quasi- uniform meshes, the assertion follows immediately from the general approximation result 
given in Theorem 14.11 below. □ 



Remark 3.1 We have only considered meshes of triangular elements on T. If the meshes 
contain also shape-regular parallelogram elements (i.e., affine images of the reference square), 
then a priori error estimates of Theorem \3.1\ remain valid only in the case of the RT-based BEM. 
This is because all the auxiliary results needed for the proof are valid in this case. This, however, 
is not true for the BDM-based BEM. In particular, the arguments in the proofs of Proposition^^ 
above and Proposition 14.11 below rely essentially on the fact that the involved polynomial spaces 
form the exact curl — div sequence (on the reference element), a property the BDM-spaces fail 
to satisfy on the reference square (see, e.g., |18j ). 

Remark 3.2 We have assumed that V is a polyhedral (closed) surface. However, all arguments 
in our proofs carry over only with minor modifications to the case of a piecewise plane orientable 
open surface V. Note that in this case there are no restrictions needed on the wave number k 
to ensure the uniqueness of the solution to the EFIE; the strongest edge singularities in (|2.5p 
have the exponents 7,? = ^ (i = 1,2); the energy space for the EFIE is H (divp,r) and 
the boundary element space X/^ consists of H(divr ,T)- conforming polynomial vector fields with 
normal components vanishing on dT (see [BJ 0] for details). 



4 General approximation result 

In this section we prove the following general /ip-approximation result for the vector field u 
given by formulas (I2.2p - (|2.8p with the singularity exponents 7? and A^ (i = 1,2) satisfying the 
minimum requirements to guarantee u € X. 

Theorem 4.1 Let the vector field u be given by (|2.2p - (|2.8p on T with 7f, 7| > and X\, A2 > 

— ^ for each edge e and every vertex v. Also, let vo € V , eo € E(vq) be such that 

min{Af + 1/2, \f + 1/2, 7l e ° , 7 2 e °} = min min {A? + 1/2, A^ + l/2, 7l e ,7 2 e } 

v£V,e£E{v) 
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with X 1 - and jf (i = 1,2) being as in (|2.5p - (|2.8p . Then for any h > and for every p > 
min {X\° , A2 , 7i° - ^,72° - \], there exists u hp G X hp such that 



lu-u^Hx < C max|/i min { fc - p > +1 / 2 p" (fc+1/2) , 



, x min{Al°+l/2,A?>+l/2 l7l e ° l72 e °} ^ 

£) }• (") 

where (3 and v are defined by (13.7P and (13.8h . respectively. 

If 1 < P < min {A^'°, Ag , 7i° — 2") 72° — 2"}' ^ en ^ or an ^ ^ > i/iere exists u fcp G X/^ such 
that 

||u - u^|| x < C ^{**HV2. (4.2) 

In order to prove this theorem one needs to find discrete vector fields belonging to X/^ and 
approximating the smooth and singular parts of u such that the approximation errors satisfy 
the upper bounds in ()4.ip and (|4.2p . 

We start with formulating the following /ip-approximation result for regular vector fields on 
r. This result will be used, in particular, to approximate the vector field u reg G X fc in (|2.2p . 



Proposition 4.1 Let P^p : X — > X/j p be the orthogonal projection with respect to the norm in 
X. If u G X fc with k > 0, then 

||u - PftpUHx < C /t min {fc.p}+l/2 p -(*+i/2) || u || xfc (4.3) 
with a positive constant C independent of h, p, and u. 

The proof is given in [4, Theorem 4.1] for the case of RT-spaces, and it carries over without 
essential modifications to the case of BDM-spaces on triangular elements (cf . Remark 13. ip . 

Now, we will study approximations of the singular part u s j ng in representation (|2.2p . By 
(|2.4p - ([2.8p and due to the arguments in Remark 12.21 (i), we conclude that u sing can be written 
as 

u S mg = curl r w + v = curl r w + (v±, v 2 ), (4.4) 

where w G H l / 2 {Y) and v G H 1 / 2 (r). 

Let us define the following discrete space (of continuous piecewise polynomials) over the 
mesh A^: 

S hp {T) := {v G C°(T); v\ Fj o Tj G V P (K), j = 1, . . . , J}. 
We will also need the following functions of h and p: 

f^Kp) := hVp-^O- + ]og(p/h))^, (4.5) 

where j = 1,2, 

ctj :=mm{Xf + l/2, 7j e °}, (4.6) 
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{ n Vo _i_ „eo i 1 ;f \Vo _ _ eo 1 
1j + S j + 2 -7i -2> f4?) 

gj° + sf otherwise, 1 ' 

otherwise, 

and the numbers 7J , AJ°, s^ , gj° (j = 1,2) are given in (|2.5p — (|2.8j) for the vertex-edge pair 
(t>o,eo) introduced in the formulation of Theorem 14.11 

In the next two lemmas we formulate approximation results for the vector fields curlr w and 
v on the right-hand side of (|4,4p . The proofs are given in Section [5] below. 

Lemma 4.1 Let w € H l l 2 (Y) be the scalar singular function in representation (|4.4p . Then for 
any h > and p > 1, there exists w hp S Shp(T) such that curlr u"^ € ~Khp and there holds 

11 1 1 ton ^\ C h{h,P) if P>ai-\, 

curl r w - curl r w p x < < , ,„ •, (4.9) 

" {ChP +1 / 2 if l<p<a x -\, 

where fi(h,p) and a\ are defined by (j4.5|) and ()4.6|) . respectively. 

Lemma 4.2 Let v £ H 1 / 2 (T) 6e the singular vector field in representation (|4.4p . Then for any 
h > and p> 1, there exists v hp £ X^p satisfying 

>, \ C f2(h,p) if p > a>2 — h, 

llv-v^Hx < < l',, 1 (4-10) 

where f2(h,p) and 02 o^e defined by (j4.5f) and (|4.6p . respectively. 

Now we are able to prove Theorem 14. 1[ 

Proof of Theorem 14.11 For the regular vector field u rcg € X fc in (|2.2p we use the orthogonal 
projection : X — > X^p with respect to the norm in X to define UrJg := P/i p u re g £ X^ p . Then 
we have by Proposition 14.11 



Keg " U&||x < C h min ^ +1 l 2 p-( fe + X /2). (4.H) 

Since the singular part of decomposition (|2.2p can be written as in ()4.4p . we use approximations 
w hp € S'^p(r) and v hp £ X^ p from the above two lemmas to define u^ g := curlr w hp + v hp £ 
X/ip. Then, applying the triangle inequality, we obtain by (|4.9p and (|4.10p 



■ I hp „ x }Cmax{/i(/i,p), / 2 (/i,j))} if p > min{Qi,Q 2 } - 5, 

lining " U si ^ ||x < < /2 . x (4.12) 

6 [(7W +i/z if 1 < p < mm{ai,a 2 } - 5. 

Now, we set u hp := UrJg + u^ g £ X/j p . Combining estimates ()4.1ip and (|4.12p . applying 
the triangle inequality, and using expressions (|4.5p and (|4.6|) for the functions fj(h,p) and the 
parameters a,-, respectively, we prove the desired estimates in (|4.ip and (|4.2p . □ 
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5 Proofs of technical lemmas 



In this section we prove Lemmas 14.11 and 14.21 
5.1 Proof of Lemma 14.11 



Recalling Remark 12.21 (ii) , we use the results of [31 |T| to find the desired piecewise polynomial 
w hp G S hp (T) such that the norm \\w — w p \\ h i/2(y) is bounded as in (14. 9p (see Theorems 5.1, 5.2 
in [3] and Theorem 4.1 in [1]). Then, recalling the fact that the the operator curlr : H l / 2 (T) — > 
Hm (r) is continuous (see |12j), we derive the estimate in (|4.9|) : 



| curlr w — curlr w hp || x = || curlr (tu — w hp )\\ l , i -> 



H7 1,a (r) 



< C\\w - w hp \\ H i/2, r) < 



C fi(h,p) if p>ai-L 



2(r) -{ChP +1 / 2 if l<p< ai -\. 
It remains to prove that curlr w hp € X^ p . In fact, it is easy to check (see [61 p. 615]) that 
A^Ti(curl r w hp \ Fj ) = curl (w hp \ T] o Tj), where curl = 



Hence, recalling that w hp \ Vj oTj G V P (K), we conclude that Mj l (curl r w hp \ Vj ) G (P p -i(K)) 2 C 
V P {K) (this is because (P p _i(iT)) 2 is a subset of both Vp DM (K) and -pJ^iT)). Moreover, 
curlr w hp G H(divr,T), because divr (curlr w hp ) = on F. Therefore, curlr w hp G X^ p , and 
the proof is finished. 

5.2 Proof of Lemma 14.21 

Let v = (vi,v 2 ) be the second term in decomposition (|4.4p of the singular vector field u s j n g. 
Again, using the results in [HI H], we find continuous piecewise polynomial approximations to 
the scalar components v±, v 2 of v: for any h > and every p > 1, there exist vl P , G Sh p (T) 
such that for i = 1, 2 there holds 

11 ftpii ,{Cf 2 (h,p) if p>a 2 -±, 

pi - v i p \\ H i/2 {r) < < , 1 (5.1) 

v ' y Ch p+i/z if 1 < p < a 2 - ^ 

with positive constants C > independent of h and p. 

Let v hp = (wf p , v 2 p )- Then v hp G H(divr,T). Furthermore, in the case of BDM-elements 
we observe that for any Tj there holds 

Mf(y hp \ v .) = J.Bj 1 ^^) o Tj G (V p {K)f = Vf M {K). 
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Therefore, v hp G X/y, in this case. Moreover, since v G H 1 / 2 (r) and v^ p G S , /i P (r), we use 
estimate (|5.ip and the continuity of the operator divr : H 1 / 2 (r) — > H~ l l 2 (T) to obtain 

||v-v^||x = ||v-v^|| HFl/2(r) + ||div r (v-v^)||^ 1/2(r) 

2 

< C||v-v h ' p || Hl/2(r) < C^2\\ Vi -t;f p || H i/2 (r) 

i=l 

< (Cf 2 (h,p) \i p>a 2 -\, 

- if i <p< a2 _ i. 1 • > 

Unfortunately, in the case of RT-elements this component-wise approximation of v does not 
work since the dimension of the RT-space on the reference triangle is smaller than the dimension 
of the BDM-space, and, in general, -MjV^lr, ) t V* T {K), so that v hp = (v^, v^ p ) <£ X hp . 
In this case we follow the procedure described in [3] (for the scalar case) . More precisely, we use 
appropriate /i-scaled cut-off functions and represent the vector field v as the sum of a singular 
vector field if with small support in the vicinity of the edges and a sufficiently smooth vector 
field rj). In the rest of this subsection we will demonstrate how this procedure works in the 
vector case. We will give a concise step-by-step outline of the procedure referring frequently to 
[3] for particular error estimates and other technical details. 

Step 1: decomposition of v. Let ho = (aia 2 )~ 1 h with <7i, o 2 from (|3.ip . Using the scaled 
cut-off functions xl( x e2//io) and x v ( r v/ho) with x 2 an d x" from (|2.5p and (|2.6p . respectively, 
one can decompose v as follows (cf. [U eqs. (5.4), (5.21), (6.4)]) 

v = (f + ip, (5.3) 

where supp<£ C LS v£ v ^eeEtv) iAe u Aj), <P G H 1 / 2 (r), ip G X m (see (|2.3p for the notation), and 
xj) vanishes in small (/i-dependent) neighbourhoods of each vertex and each edge of V. Here, A e 
is the union of elements at one edge e, i.e., A e := U{fj; Tj n e ^ 0} (note that the endpoints of 
e are not included in e), A v is the union of elements at a vertex v, i.e., A v := U{Fj; v G Tj}, m 
is sufficiently large and depends on the parameters t 2 and m 2 specified for the singularities 
and u!f , respectively. 

Step 2: approximation of <p for p = 1. If p = 1, then one can approximate <p by zero. 
One has <p hp = G X^p, and, recalling the continuity of the operator divr, we derive 

\\<P ~ <Php\\x- = llv»ll H -i/2 (r) + " divr vWh-i/ht) < CH^IIni/afr). (5.4) 

We will obtain /i-estimates for the norms of ip in H s (r) with s = and s = ^ + e (for sufficiently 
small e > 0) by using the fact that (p has a small support. First, we apply Lemma 3.1 of [2] 
and Lemma 3.5 of [3j to localise the norm ||v?|| H i/2+ e (r) to the faces C T and to the elements 
Tj C T? , respectively. Then, we use scaling on each element Tj C supp<£> (cf. [33 Lemma 3.1, 
eqs. (5.6), (5.11), (5.23), (6.7)]). As a result we have 

[|¥>||H»(r) - C /j H^llH^r/) 
/ : r/cr 
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^ C E E f'* ^llvlli^cr^ + MWi)) 

< Ch 2 ^ +l / 2 ' s \l + \og{l/h)) 2 ^ for s G {0,1/2 + e}. 

Here, c*2 and fa are denned by (|4.6p and (|4.7p . respectively. Hence, using the interpolation 
between H°(r) and H 1 / 2 "^^), we obtain by (J53J) 

\\<P ~ ¥» hp ||x < C7^(l + log(l//i))&. (5.5) 

Step 3: approximation of <£> for p > 2. We approximate each component tpi of by a 
piecewise polynomial ip^ G Sh, P ~i(T) (i = 1,2). Here we can use the results of [3J [1] for each 
type of singularity: there exist ip\ v G 5/ l)P _i(r) such that for i = 1,2 there holds (cf. [3j eqs. 
(5.12), (5.13), (5.24), (6.9)]) 

\\m - ^ P |l^/ 2{ r) < Ch^ (p - l)" 2 ^ (1 + log 2=1)^ < Ch^p' 2 "' (1 + log I f\ (5.6) 

where a 2 and /?2 are the same as in (|5.5p . 

Setting ip hp = ((fi P , ^2 P ) we observe that <p' ip € H(divr,T), and for any element I\ one has 

MjHfPlr,) = J 3 Bj\^\ Tj )oT 3 G (V p -i(K)f c 7>f (K). 

Hence, <^ hp G X hp . Moreover, since tp G H x / 2 (r) and $ p G S hiP _i(r) for i = 1,2, we estimate 
by analogy with (|5.4p and using (|5.6j) : 

llv - Vfepllx < C\\cp - ¥> ft p|| H i/2(r) 

2 

< C^||^-^ p ||^/ 2 (r) <C^>- 2Q2 (l + log|)^. (5.7) 
i=i 

Comparing (|5.5p and (|5.7j) we observe that one can use estimate (15 .7p also in the case p = 1. 

Step 4: approximation of ^. Recalling that ip G X m we apply Proposition 14.11 there 
exists ip p G X/ip such that 

H^-^llx < C/i min{A!,p}+1/2 p- (fc+1/2) HVIIx*, 0<A;<m. (5.8) 

We also recall that ip vanishes in the /i-neighbourhood of all edges and vertices of T. Therefore, 
having explicit expressions of the singular vector field v (and thus, of ip) on each face C V, 
one can derive upper bounds for the norms 1 1 "0 1 1 n fe (r/ ) an d ||div r /t/>||#fc( r /) in terms of the mesh 
size h, the parameter k, and the singularity exponents. For the norm of ip this can be done 
component- wise using the same calculations as in [3] (see, e.g., inequalities (5.15) and (6.10) 
therein) : 

\MuHrf) < Ch a > +1 ' 2 - k {\og{\/h)f* +!> \ a 2 + \ < k < m, (5.9) 
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where oli and are defined by (|4.6p and (|4.7p . respectively, = 5 if ^ = «2 + 5 and ^2 = 
otherwise. 

To estimate the norm ||div r /^/>||#fc/ r A we observe that the operator div r / reduces all singu- 
larity exponents by one (while preserving the structure of the corresponding singularity). Then, 
using similar calculations as indicated above, we obtain 

||div r /V||/j*(r/) <C/ t Q2 - 1 / 2 -^(log(l// i ))^+-2 ) a2 _i< fc < m) (5.10) 

where 02, /?2 are the same as in (|5.9|) . whereas T>2 = \ if k = 012 — \ and V2 = otherwise. 
By (O)) - (pn01) we conclude that 

\\rl> " V^llx < C h min{k ' p}+a2 - k p- ik+1/2) (log(l//i))^ 2+p2 , a 2 - \ < k < m (5.11) 

with the same 02, P2, an d ^2 as in (|5,10p . 

Let p > 2a2 — \- Since m is large enough, we can select an integer k satisfying 

2a2 — \ < k < min{m,p}. 

Then mm {k,p} = k, and p~( fe+1 / 2 ) < p~ 2a2 . 

If a.2 — \ < p < 2a2 — 5 (i ,e- > P * s bounded), we choose an integer k £ (02 — 2>Pb an< ^ 
if p = «2 — then we take k = p = «2 — \. In both these cases min{/c,p} = k, and 

p -(fc+l/2) < Cp~ 2a2 . 

Thus, for any p > a.2 — 5, selecting as indicated above we find by (I5.11|) 

||^ - *J> hp \\x < Ch a2 p~ 2a2 (log(l/h)f 2+ » 2 (5.12) 

with Q2, /?2 , and being defined by (|4.6p . (|4.7p . and (|4.8p . respectively. 

Step 5: approximation of v = ip + -0. Let us define v hp := cp hp + ip hp £ X^p, where ip hp 
and are approximations constructed above (see Steps 2-4). Then combining estimates (|5.7p 
and (|5.12p and using the triangle inequality we prove (|4.10p in the case p > 02 — \ ■ 

It remains to consider the case 1 < p < a.2 — \ - In this case one does not need decomposition 
(|5.3p . Observe that for each face one has div r / v £ H k (F^) with 1 < k < a — ^. Therefore, 
v € X fc with 1 < k < a — i, and applying Proposition 14. II we find v hp € X/^ satisfying 

||v - v hp ||x < Ch min ^ +l / 2 ||v|| xfc . 

Hence, selecting /c € [p, a — ^) we arrive at the desired upper bound in (I4.10p . and the proof is 
finished. 
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